Disorder-induced magnetooscillations in bilayer graphene at high bias 
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Energy spectrum of biased bilayer graphene near the bottom has a "Mexican-hat" -like shape. For 
the Fermi level within the Mexican hat we predict that, apart from conventional magnetooscillations 
which vanish with temperature, there are additional magnetooscillations which are weakly sensitive 
to temperature. These oscillations are also insensitive to a long-range disorder. Their period in 
magnetic field scales with bias, V, as V^. The origin of these oscillations is the disorder-induced 
scattering between electron- like and hole- like Fermi-surfaces, specific for Mexican hat. 

PACS numbers: 73.21.Ac,73.20.At,73.43.Qt,71.55.Jv 



Introduction. Numerous experimental studies of elec- 
tronic properties of bilayer graphene were recently re- 
ported in the literature |ll4l5|. From prospective of po- 
tential applications, the appeal of bilayer graphene is that 
a gap can be opened and tuned by the gate voltage |6l4l5j . 
Using the dual (top and back) gated structures allows 
to control both the gap and the carrier density indepen- 
dently |lll4l5l |. For these structures, opening of a gap was 
demonstrated in temperature [ll| and bias [ll|, |l2[ de- 
pendence of resistivity, in bias dependence of capacitance 
14 . Il5| , as well as in strong- filed magnetotransport tl3| . 
Measurements in magnetic field reported in the literature 
focused either on weak-field {B ~ 0.1 T) domain in or- 
der to reveal weak localization [2j and universal conduc- 
tance fluctuations [lO[, or quantizing {B ^ 10 T) fields 
SSEIillil- At intermediate fields, B - 1 T, trans- 
port and capacitance are determined by electron states 
near the band-edge, Fig.[T] As follows from tight-binding 
calculation by McCann and Fal'ko [l6|, the spectrum 
near the band-edge has a form of "Mexican hat" with 
minimum at 
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where v = 8 x 10^ cm/s is the band velocity, t is the 
interlayer hoping and V is the bias. The minimum has a 
depth. 



FIG. 1: Energy spectrum of a biased bilayer graphene near 
the bottom of conduction band. Two Fermi circles, electron- 
like and hole-like, have radii pe and Ph, respectively. Black 
dots are electron-like Landau levels. White dots are hole-like 
Landau levels. Red arrows illustrate scattering process giving 
rise to additional oscillations. Green arrows illustrate that 
energies of electron-like Landau levels, Eq. ((8)1, grow with 
number, while energies of hole-like Landau levels fall off with 
number. 



Ejn, the behavior of magnetocapacitance, SC{B)/C^, and 
magnetoresistance, SR{B)/R, exhibits oscillations, which 
are additional to conventional magnetooscillations re- 
lated to the Landau levels. They behave as 
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The "capacity" of the minimum for t = 0.4 eV and 
V = 100 mV is n = pl/{TTh^) = 0.54 x lO^^ cm'^. This 



density is comparable to the densities in experiment IJ] , 
but the gap in this experiment was as small as 26 meV. In 
experiment 15| the gap was wider, w 100 meV, but the 
density was higher, ~ 1.5 x 10^^ cm~^. In both experi- 
ments, the Fermi energy exceeded Em- Upon increasing 
the bias from V = 100 mV to T^ = 250 mV, as in exper- 

This 



12l |. the capacity increases by a factor 4.8 



iment 

suggests that in the situation Fig. [1] when all electrons 

reside in is feasible. 

In the present paper we demonstrate that, in the 
regime Fig. [l] when the Fermi energy is smaller than 



where A = y^hc/{eB) is the magnetic length and po is 
the minimum position of the spectrum, given by Eq. ([T|) . 
The remarkable feature of the oscillations Eq. ([3]) is that 
their period does not contain the Fermi energy, Sp. This 
suggests that they are not smeared upon increasing tem- 
perature, while conventional oscillations are suppressed 
as exp{—27r^T/hujc)- For the same reason, additional os- 
cillations Eq. ([3]), unlike conventional magnetooscilla- 
tions, are insensitive to the random density variations 
caused by long-range disorder. On the other hand, oscil- 
lations Eq. ([3]) are disorder-induced, since it is scattering 
by short-range disorder which causes the term Eq. ([3|) in 
the response functions. 



The origin of additional oscillations lies in the fact that, 
for Sp < Sm, the Fermi surface consists of two circles 
with radii, Peisp) and Phisp), Fig. [1] Magnetooscilla- 
tions corresponding to these circles are cos(7rpgA^) and 
cos(7rp^A^), respectively. Due to disorder scattering be- 
tween the two Fermi circles, magnetocapacitance and 
magnetoresistance will also contain a product. 
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cos(7rp2A2) cos(7rp^,A2) = - cos(7r [pi - pi] X 
+ -cos{tt [pI+pI]X''). 
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Our prime observation is that Ep drops out of the second 
term. Indeed, the spectrum of bilayer graphene has a 
form 
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so that 
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We see that relation pi + pf, = 2pQ holds for arbitrary 
£p < Em- In the reminder of the paper we give a deriva- 
tion of Eq. ([3|) with a prefactor and briefly discuss the 
limit of strong magnetic fields, where the Hall quantiza- 
tion becomes important. 

Density of states. We start from zero magnetic field. 
Spectrum Eq. ([5]) is a result of diagonalization of the 
2x2 matrix Hamiltonian 
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It is convenient to find the Landau level spectrum of Eq. 
d?]) using the gauge, A = (0, Bx), in which eigenfunctions 
of the conventional quadratic spectrum are e*P»^(/>„(a;), 
with (f>nix) being the eigenfunctions of ID harmonic os- 
cillator. In this basis, the spectrum for n > 2 is 
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with Q.C ~ 2hv /(tX ), while eigenfunctions are given by 
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where 



dn = 



Mlc^n{n — 1) 



(10) 



Spacings between the electron-like and hole-like Landau 
levels at the Fermi energy are the same. 
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In the presence of disorder, U(r), imaginary part of 
the self-energy, 'E,n{E), is determined by level numbers, 
n, for which e„ is close to E. At the same time, the n- 
dependence of E„(i?) is weak. In our case, however, the 
equation Sn — E has two different solutions: 2ne ~ P^X"^ 
and 2n/i ~ Ph^^ ■ Correspondingly, one should consider 
two different self-energies, I]„^ = Se and I]„^ = S/j. As 
demonstrated in [19|, short-range disorder (with corre- 
lation length smaller than A) insures the applicability of 
the self-consistent Born approximation (SCBA), which in 
our case becomes a system 
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where subscript (n, e) or (n, h) indicates that the sum- 
mation is performed for n close to Ue or n/j, respectively. 
Coefficients 
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do not contain numbers of Landau levels, since they are 
taken at n = ne and n = rih, correspondingly. Explicit 
expressions for dn^ and d^^^ as a function of the Fermi 
energy are the following 
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where em is given by Eq. ^, and Pe, Ph are given 
by Eq. ©• Coefficients ath and ate in Eq. P^ are 
given by Eq. (J14p with replacement, e ^ h. Coeffi- 
cient r2 describes the strength of the disorder potential, 
r2 = /d2r([/(o)[/(r))/(27rA2). Second terms in Eqs. 
([T2l) . (fT3|) describe contributions from disorder-induced 
mixing of electron-like and hole-like Landau levels. 

Applying the Poisson summation formula to the sums 
in Eqs. ([T^. p^ . we get the following equations for the 
imaginary parts ImE = S", 
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where Tg and t^ are the full scattering times from the 
states Pe and ph in a zero magnetic field, 
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We note that magnetic-field dependence drops out from 
the ratio V"^ / {fiWc)- Concerning the energy dependence of 
Te and r/i, it follows from Eq. ((TT|) that near the bottom 
of Mexican hat, Ep <C £m, we have Te, r/j ~ v^Ef, which 
refiects the ID character of the bear density of states [20j. 
Iterating Eqs. (TTSl) . (flT)) . we obtain a contribution to the 
density of states of the form 



5g{B) = Go cxp 



UcTe UJcTh 
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which coincides with additional oscillations Eq. Q 
stated in the Introduction. Prefactor Gq is given by 
Go = inT'^aeh/ [(^c)^A^] . Magnetic field dependence 
of Go is oc 1/B. Energy dependence of Go is plotted in 
Fig. [2] We see that Go diverges in the limit e^ ^ 0. 
This divergence is also due to the ID character of the 
density of states near the bottom of the Mexican hat. 
As the Fermi level approaches the top of the Mexican 
hat, the hole contributions in Eqs. ([TS]). p7)) . and result- 
ing additional oscillations, disappear. At the same time 
the prefactor Go remains finite in the limit Ep — > £,„. 
Such an abrupt behavior of additional oscillations is a 
consequence of the fact that the bare density of states 
experiences a jump at Sp = Em- 
Conductivity To trace the emergence of the product 
cos(7rpgA^) cos(7rp^A^) in the conductivity, a, it is suf- 
ficient to set Eg = ft/re and EJJ — h/rh in the exponents 
in Eqs. (HH), (HZ]). The SCBA expression for a in the 
case of bilayer graphene is the sum of electron-like and 
hole-like contributions 

a{E)^a,iE)+ahiE), (20) 

where ae can be presented in the form of a sum, 
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and ah is given by Eq. (|2ip with replacement of subindex 
e by ft. Matrix element {vx)n.n+i taken between the 
states Eq. © is 
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FIG. 2: Prefactor in "magnetocapacitance" 5g{B), Eq. (|19|l . 
in the units of 47rr^A'' [2pot/{h^vV)y\ is plotted from Eqs. 
(Ilip . (|14p . and (|15|l . versus dimensionless ratio ep/em for 
three values of the ratio, V/t = 0.5, V/t = 0.7, and V/t = 1. 



Now we notice that a term proportional to the product 
of two cosines Eq. ^ follows from (S") in the numer- 
ator of Eq. (HJ). Indeed, plugging Eq. p^ into the 
numerator of Eq. (|2ip and replacing (S") in denomina- 
tor with its leading value ft^/r^, we arrive at additional 
oscillations 
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where we have introduced the Drude conductivity 
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A more accurate form of Sae/u'^ can be found by applying 
the Poisson summation formula to Eq. (j2ip . which will 
include corrections to self-energies in the denominator. 
This transforms Eq. (1^51) into 
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which differs from Eq. (I23P in strong, WcTg > 1, magnetic 
fields. For the contribution ah, relations Eq. ((25)) holds 
with replacement of subindexes e ^ ft. 
Concluding remarks, (i) Up to now, observing the Mex- 
ican hat structure of the spectrum in experiments on 
bilayer graphene was limited by relatively low mobility, 
fi - 2000 - 3000 cm2/Vs. This corresponds to the en- 
ergy smearing, ft/r, of about 7—10 meV. In particu- 
lar, revealing Landau quantization in magnetotransport 
[ll lD, lid ll3| , ac [J] , and magnetocapacitance IJ, ll5[ ex- 
periments required strong magnetic field, i? ^ 10 T. Typ- 
ical cyclotron quantum for such fields is hcoc ^ 20 meV, 



i.e., it is bigger than £„ « 11 meV for V = 0.2 V . On 
the other hand, inhomogeneity of local electron density 
[lj| was a significant factor in smearing of magnetooscil- 
lations. In this regard, additional oscillations Eq. ([3]), 
being insensitive to this inhomogeneity, might be ob- 
servable even when conventional magnetooscillations are 
completely washed out. For B = IT and same V — 0.2 V 
condition Em > fy^c is satisfied; for this V , the product 
PgA^ in the argument of Eq. ([3]) is 38. For such fields, 
conventional oscillations are suppressed at temperatures 
as low as r = 3 K, while additional oscillations Eq. ^ 
remain unaffected. Unlike conventional magnetooscilla- 
tions [13| , they are also insensitive to the lifting of valley 
degeneracy. 

[ii). Our calculation was based on the spectrum 
Eq. ([5|); this spectrum is obtained from 2x2 Hamilto- 
nian Eq. ([7]) . Analysis of more general 4x4 Hamiltonian 
16|,ll7| suggests that the gap can exceed t while the prop- 
erty, pI + pI = 2pq persists. 

[iii). Relevant densities for the additional oscilla- 
tions are ^ 10^^ cm^^. Such densities are high enough 
for electron-electron interaction- induce d spectrum renor- 
malization to be insignificant 



22|-|24{. On the other 



hand, interactions can scatter electrons between electron- 
like and hole-like Fermi surfaces. They also induce in- 
elastic lifetime ~ Sp/T'^. This leads to effective sup- 
pression of additional oscillations at temperatures above 
T - y/eju^c ~ 50 K . 

[iv). To establish a relation between oscillation Eq. ([3]) 
and magnetointersubband oscillations in a quantum well 
with two subbands [2l| , let us turn to the product Eq. (j4]) 
of the oscillating part of the density of states. Magne- 
tointersubband oscillations of Ref. [2l| follow from the 
similar product for different subbands. However, they 
emerge from the term, cos[7r(pg— p^)A^], of Eq. (jl]), while 
oscillations Eq. ([3]) come from the term cos[7r(pg -\-p\)}?] 
of Eq. (j4|). Independence of this term of e^ is specific for 
bilayer graphene. 

[v). In closing, we discuss qualitatively the limit of 
quantizing magnetic fields. When the Fermi level lies 
within the Mexican hat, classical trajectories correspond- 
ing to electron-like and hole-like states are Larmour cir- 
cles with opposite direction of rotation. Indeed, the equa- 
tion of motion in momentum space, p = -4^ x B, can 
be presented as 
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With energy and absolute value of momentum conserved 
by Eq. (1^5)) . the only difference between electron-like and 
hole-like motions comes from the factor, (p^ ~Po)- Since 
{pI — Po) — ~{Ph — Pq), clockwise rotation of electrons- 
like states and anti-clockwise rotation of hole-like states 
have the same frequency, in agreement with Eq. ()11|) . At 
the same time, the radii and velocities of their Larmour 
motions are related as Pe/ph- 



Opposite directions of rotation for electron-like and 
hole-like states translate into the opposite sings of drift 
velocities for corresponding edge states. 
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This, in turn, means that dispersion laws for electron- 
like and hole-like states intersect each other. Previously, 
Refs. 25|, l26| pointed out that opposite dispersion of 



the edge states from the same Landau level can arise 
from the valley splitting. Combined with the Zeeman 
splitting, this leads to intersecting edge dispersions for 
opposite spin directions 25|. We note that in bilayer 
graphene with the Fermi level within the Mexican hat, 
crossing of the edge dispersions from different Landau 
levels occurs naturally with the valley degeneracy pre- 
served. For the Fermi level located at the intersection 
of electron- and hole-dispersion curves, interactions can 
result in non-chiral Luttinger liquid at the edge. This sit- 
uation is similar to the quantum Hall line junction con- 



sidered in Refs. [27|,l28|. Unlike Refs. [271, l28i , where the 



disorder results in resonant-tunelling states between the 
edges separated by a tunnel barrier, in our case disorder 
will smear the corresponding Luttinger-liquid anomalies. 
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